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I . INTRODUCTION 


This report describes an anlaysis of the small scale wind velocity 
perturbations in vertical wind profiles at Cape Kennedy, Florida. The 
overall objective of the analysis is the derivation of 'new information 
that is useful for simulations of Space Shuttle ascent thru the perturbed 
atmosphere. Previous representations of wind velocity perturbations 
include power spectra (Ref. 1, 2, 3), exceedance statistics (Ref. 3) and 
idealized singularities and quasi-sinusoidal wave patterns (Ref. 2). 

Each of these representations by themselves is an important aspect of 
perturbation profiles. Power spectra are suitable for representation of 
the random component of the perturbations, Perturbations having a known 
power spectrum are simulated by application of an appropriate filter to 
a noise generator. These simulated perturbations do not contain the 
singularities and quasi -sinusonidal variations that are observed in wind 
profiles, Prior to this study, the available statistical data did not 
permit specification of various aspects of idealized singularities and 
wavelike perturbations with a reasonable degree of confidence. The 
information developed as a result of the analysis described in Section III 
of this report is suitable for the further development of idealized 
models. 


The term perturbation is used here instead of the more common 
term, gust. According to the conventional approach, a gust profile is 
calculated by applying a high-pass digital filter to a Jimsphere profile; 
all the speeds in the filtered profile are defined as gusts. In this 
study the high-pass filtered profile is defined as a residual profile 
and the maximum residual 1n the vicinity of a specified reference height 
is defined as the gust. Gusts defined in this manner represent the 
perturbation peaks. A detailed discussion of the calculation of residual 
profiles and gusts is given in Section II. The meteorological coordinate 
system, the data sample and Jimsphere profiles are also described in 
Section II. Recommendations and conclusions are presented in Section IV. 


II. TECHNICAL BACKGROUND 


A. Meteorological Coordinate System . 

The basic winds aloft data are recorded in terms of wind direc- 
tion, 6 and magnitude, W. The wind vector is expressed in the standard 
meteorological coordinate system in which the direction from which the 
wind is blowing is measured in degrees clockwise from true north. The 
zonal component, u, of the wind vector is positive for a west (west to 
east) wind (6=270°) and negative for an east (east to west) wind (6»90O); 
the meridional component, v, is positive for a south (south to north) 
wind (0=180°) and negative for a north (north to south) wind (9*0°). 

The components are calculated from 0 and W according to; 

u = -W sin 9, 0 < 6 < 360° (1) 

V = -W.cos 6, (2) 

The relation hntmen 6 defined above and the angle defined in 
standard mathematical polar form is: 

0 = 270 (3) 


B. Data Sample . 

The data consist of 1800 Jimsphere profiles (160 per month) 
from Cape Kennedy, Florida (Ref. 4). The data were obtained under a 
Space Shuttle Level II directive that specifies the demonstration of 
vehicle design validity using 150 Jimsphere wind profiles representative 
of each month. Two months (April, February) were chosen for analysis 
in this study. April data were used to develop and refine the analysis 
procedure which could be applied efficiently to other months when 
required. April was also of interest because it coincided with the 
planned(l) Orbital Flight Test Mission. The February dcta were chosen 
because they are representative of the winter season at Cape Kennedy. 

The number of soundings for each month for each year of the sampling 
period is illustrated in Figure 1. 

C. Jimsphere Data . 

Jimsphere wind profile data are obtained by precision radar 
(FPS-16) tracking of an ascending (4-5 m/s) 2-meter diameter super- 
pressure balloon with a roughened surface. The balloon positions, deter- 
mined every 0.1 second, are smoothed to provide mean positions at each 
25-meter interval of ascent (Ref. 3). Differences in position between 

^^^Rescheduled 
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Figure 1. Distribution of February and April Jimsphere Soundings 
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alternate 25-meter levels then indicate the mean wind for the corresponding 
50-meter, layer, and are reported as the wind at the 25-meter level in the 
middle of the 50-meter layer. Thus, the basic data analyzed here are 
wind speeds and directions for 50-meter layers, overlapping by 25 meters. 
Even when the overlaps are eliminated, winds for successive 50-meter layers 
are not independent, because they are based on the smoothed balloon 
position at the common boundary. Only when at least 25 meters intervenes 
between two layers (i.e., winds reported for levels at least 75 meters 
apart) can two winds be considered independent observations (Ref. 3). 

Expressions for the amplitude response, G(X), of the Jimsphere 
system to wind perturbation wavelengths that are small relative to the 
length of the wind profile have been derived by Luers and Engler (Ref. 5), 


G(X) * 



and by DeMandel and Krivo (Ref. 6). 


(4a) 



(4b) 


where, S = smoothing interval = 75m 
X = wavelength(m) 
w = Jimsphere balloon ascent rate 

The amplitude response is equivalent to the ratio A(X)/A*(X); where A*(X) 
is the true amplitude and A(x) is the amplitude measured with the Jimsphere 
system. As illustrated in Figure 2, the Jimsphere system does not 
measure wavelengths less than 50m; for X=90m the measured amplitude is 
one-half the true amplitude. 


The Jimsphere data represent an order of magnitude improvement in 
resolution when compared to conventional Rawinsonde wind profiles. The 
perturbations observable with the Jimsphere system permit a more realistic 
assessment of the structural response of spacecraft and launch vehicles 
to small scale wind gust and shear. 
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Figure 2. Amplitude Response of the Jimsphere System 


D. 


Filter Functions. 


Certain characteristics of winds aloft that are important in 
the planning and mission operations of space vehicles launched or retrieved 
at a particular location can be singled out for detailed study byapplica- 
tion of specially designed digital filters to Oimsphere wind profiles. 

The design of the digital filters is based on the Martin-Graham 
cosine rolloff model described by Demandel and Krivo (Ref. 7). A set of 
numerical smoothing weights is calculated for a low-pass filter from the 
equation 


h(nT) 


sin(2iTf^nT) + sin(2iTf nT) 


2irnT 


L 1 - - fc>'J 


(5) 


Where the filter design parameters are 

T = altitude interval of wind profile data 

n = weight index (-N, -N+1,.. . ,-1,0,1,. .. ,N-1,N) 

N = (NW-l)/2 

NW = number of weights 

f = cutoff frequency = the highest frequency with associated 
amplitude passed with unity gain 

f^ = termination frequency = the lowest frequency with 
associated amplitude passed with zero gain. 

The center weight (n = 0) is given by: 


When the weights, h , have been determined, they are normalized by 
applying the constraint 

N 

( 6 ) 

n=-N 

Only (N + 1) weights are calculated since hn = h_n. Since the filter 
function is symmetrical, no phase shift is produced. 

The use of digital smoothing weights results in the loss of the 
first and last N data points of the original profile. Thus the filtered 
wind profile has an altitude range that is reduced by 2NT compared to 
the original profile. 
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The effective response of the filter, given the design paramete*'? 
listed under equation (5) is 


N 

S(i^) = ^ + 2 2 h^ cos(2TrfnT) , (7) 

n=l 


As the number of weights (NW) is increased, the response of the filter 
improves. However, computation time increases as does the number of 
points lost (the first and last N data points). In this study NW was 
chosen to minimize data loss while maintaining a reasonably accurate 
filter response. 


E. Filter Application . 

Jimsphere wind profiles in component form (zonal and meridional) 
were decomposed into five data bases by the filtering process illustrated 
in Figure 3. The filter weights, h^, and gain functions, G(f) calculated 
from equations 5 and 7 for low-pass filters I and II are listed in Tables 
1 and 2, respectively. The power response of the filters, G2(f), illustrated 
in Figure 4 is a measure of the effect on the variance of the perturbations 
in the profile. For example, low-pass Filter I reduces the variance to 
39.5 percent of its original value at a wavelength of 500 meters. 

The effect of the filtering process on the zonal and meridional 
components of a particular profile is illustrated in Figures 5 and 6. 

The terms (steady state, residual I, wind bias, residual II and total 
residual) used to describe the data bases generated by the two-stage 
filtering process are defined below. 

Steady state profiles are appropriate for analysis of winds aloft 
OP the synoptic scale; these profiles represent average conditions over 
spatial scales of a few hundred kilometers in the horizontal and a few 
hundred meters in the vertical and temporal scales of a few hours. 

The first residual profile (residual I) is calculated by subtraction 
of the steady state profile from the Jimsphere profile. This is analogous 
to the application of a high-pass filter having an amplitude response 
function that is the complement of the low-pass filter used to calculate 
the steady state wind profiles; i.e., 

G(f) high-pass = 1 - G(f) low-pass (8) 


The first residual profile contains the small wavelength perturbations 
in the Jimsphere profile that are not found in the steady state profile. 
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Figure 3. Digital Filtering of Oimsphere Profiles 






TABLE 1, FILTER WEIGHTS, hi, AND GAIN FUNCTION, ,G(f) OF LOW-PASS FILTER I 
FOR T»25 tn, N»20, fc“. 00034 nr^ AND 00435 fn“l 




Filter 




Altitude 


Weight 

X(m) 

G(f) 

2 (1) 

h- 

.116360 

10,000 

.0001 

.9994 

+ 25 m 

0 

^+1,-1 

.112682 

5,000 

3,333 

.0002 

.0003 

.9977 

.9947 

Zq jL 50 m 

M 

^+2,-2 

It 

.102183 

.086370 

2,500 

2,000 

1,667 

.0004 

.0005 

.0006 

.9904 

.9844 

.9768 

II 

II 

.067415 

1,429 

.0007 

.9672 

II 

II 

.047750 

1,250 

.0008 

.9555 



1,111 

.0009 

.9414 

II 

M 

.029618 

1,000 

.001 

.9249 

!l 

M 

.014711 

500 

333 

.002 

.003 

.6299 

.2512 

11 

M 

.003949 

250 

.004 

.0219 

M 


-.002561 

200 

167 

.005 

.006 

-.0004 

-.0006 

U 

11 

-.005395 

143 

.007 

.0008 

II 

II 

-.005565 

125 

111 

.008 

.009 

-.0009 

.0009 

U 

II 

n 

II 

11 

M 

-.004229 

-.002423 

-.000884 

100 

.010 

-.0009 

II 

II 

II 

II 

II 

II 

.000021 

.000259 

.000022 




II 

II 

II 

-.000406 

-.000771 




Zq + 500 

^+20,-20 

-.000926 





(1) Zq = 500,525...19,500 m 
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TABLE 2. FILTER WEIGHTS, h. , AND GAIN FUNCTION, G(f) OF LOW- PASS FILTER II 
FOR T=250 m, N=5, f^=.00004^m“l AND f^=. 00080 m-1 


Altitude 


Filter 

Weight 


ho 

.20333 

1 + 250 m 
0 — 

^+1,-1 

. 18260 

Zq + 500 m 

h+2,-2 

.13008 

II 

II 

.06865 

II 

II 

.02065 

Zq + 1250 m 

^5,-5 

-.00365 


X(m) f(m~MxiO~^ Gff) 


20,000 

0.5 

.9904 

10,000 

1 

.9620 

6,667 

1.5 

.9161 

5,000 

2 

.8546 

3,333 

3 

.6955 

2,500 

4 

.5106 

2,000 

5 

.3283 

1,667 

6 

.1736 

1,429 

7 

.0627 

1,250 

8 

.0001 

’ ,111 

9 

-.0210 

1,000 

10 

-.0155 


i. 

I 





1750,1775,.. .18, 250m 


k 
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Figure 4. Power Response of Digital Filters (I and II) and Associated Band-Pass of 
Filtered Profiles (Steady State, Residual I, Wind Bias and Residual II) 
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Wind bias profiles emphasize the predominant large scale change 
in the wind from the surface to 20 km. Typically at Cape Kennedy, this 
large scale change is dominated by the maximum vector wind magnitude which 
occurs near 12 km. Steady state wind profiles clearly contain the larg-s 
wavelength perturbations needed for wind biasing, but they also contain 
small wavelength perturbations that cannot be used. Wind bias profiles 
are calculated by application of a low-pass filter to the steady state 
profiles. The calculated wind profiles contain perturbation wavelengths 
similar to that found in a wind bias model based on the monthly mean wind. 
However, large differences of wind magnitude can occur between a particular 
wind bias profile and the monthly mean profile; this is illustrated in 
Figure 7. A filtered wind profile that is representative of the wind 
conditions associated with a particular launch is the most desirable 
basis for wind biasing of launch vehicles. The monthly mean profile is 

almost never representative of launch conditions. The filter function 
used to calculate the wind bias profiles could be used in future simula- 
tion studies or launch operations. The ultimate choice of an appropriate 
filter will be based on engineering considerations beyond the scope of 
this investigation. 

The second residual profile (Residual II) contains. the small scale 
perturbations in steady state wind profiles that are not required for wind 
biasing. 

The total residual profile (Residual I + Residual II = TRP) contains 
all the perturbations that are of interest for analysis of vehicle response. 

The wavelength limits of the various types of wind profile data are 
listed in Table 3; with exceptions noted in the table, the indicated 
wavelength corresponds to the wavelength at which the amplitude response 
of the appropriate filter (high-pass or low-pass) is equal to .50. 


Table 3, Wavelength Limits for Jimsphere and 
Various Filtered Wind Profiles 


Wavelength Range (m) 


Jimsphere 

90^1) 

-X <2) 

Steady State 

420 

n,ax(2) 

'^max 

Residual I (RP I) 

90 

- 420 

Wind Bias 

2,470 

- X (2) 
max 

Residual II (RP II) 

420 

2,470 

Total Residual (TRP) 

90 

- 2,470 


^^^Valve of X for G(X)=0.5 in equation (4) 

( 2 ) 

\ ^max equivalent for Jimsphere, steady state and wind bias 
profues. The exact value of X^ax could be the subject of a detailed 

- • 


discussion which would not be pertinent to this report* If these 
profiles are subjected to a spectral analysis reliable estimates 
of the spectral energy could only be obtained at wavelengths con- 
siderably smaller than the length of the profile. However, a simple 
inspection of winds aloft statistical summaries and individual wind 
profiles leads to the conclusion that a wavelength somewhat larger 
than the length of a Jimsphere profile is detectable. 

F, Definition of Gust 


Let u' represent the zonal wind component at a specified 
r^>ference altitude, Hq, in a residual profile. The zonal gust is defined as 
the maximum value of u' in the vicinity of altitude Hq with like sign to 
u' at h'o- The wavelength, X, of the gust is defined as twice the altitude 
difference between the first zero crossing of u' detected by scanning 
upward and downward from the maximum value of u'. The altitudes of the 
zero crossings, H? and Hi, are calculated by linear interpolation according 
to 


^^2 = Hj -1 


25 


^j " ‘'j-1 


^j-1 


( 9 ) 


w. 

"k+l 


25 

“k+l ■ “k 


•k+l 




where H2 = altitude of the first zero crossing for the upward scan 

u'._, = last value of u' with like sign to sign of u' at H^ when 
scanning upward(l) 

u'. = first value of u' with sign opposite to sign of u' at Hq 

^ when scanning upward 


H.- , = altitude of ui , 

j-1 j-1 

H^ = altitude of the first zero crossing for the downward scu n 

^k+l ' value of u' with like sign to sign of u' at Hq when 

^ ^ scanning downward 

UjJ, = first value of u' with sign opposite to sign of u' at Hq 
when scanning downward 

H|^^j = altitude of Uj|,^j 

Substitution of v' for u' in the definition above yields the 
definition of the meridional gust component. 
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he indices j and k increase upward 


to 








It follows from the definition of gust adopted for this study 
that a particular gust component (zonal or meridional) 1s Independent of 
the other component. In most cases the altitude of the maximum speed 
associated with the gust and the wavelength of the gust will differ for 
the zonal and meridional components. The most realistic bivariate statistics 
of gusts and wind perturbations are obtained from the following combinations: 

• Zonal component gust, u'niax» and associated meridional 
component perturbation, v^, at the same altitude. 

• Meridional component gust, v'^a^* and associated zonal 
component perturbation, u', at the same altitude. The 
first combination is used in this study. 

Thus, for any particular profile reference height we are consider- 
ing, a wind vector which is composed of a gust or maximum of the zonal 
component near the reference height and an associated "perturbation" of 
the meridional component which is not necessarily the largest perturba- 
tion near the reference height. 

The remaining alternative is the assumption that the zonal and 
meridional gusts (u'^gx and v'„,ax) occur at the same altitude; this yields 
an unrealistically large value of the vector modulus; i.e.. 


V^^'max)^ (''‘max) ^ V^^'max^^ 




( 11 ) 


G. Gust length 

In connection with studies of gust gradient and wind shear, 
it is appropriate to define an altitude interval associated with a gust. 
This altitude interval defined as the gust length, L, is calculated by 
taking the altitude difference of the zero c^^ossings on either side 
of the gust; i.e., 

L = Hg > (12) 


III. ANALYSIS 

A. Maximum Vector Wind Modulus 


The maximum vector wind modulus, R[r,ax» '’s the largest wind 
speed, R, that can be found in a particular profile. 


(13) 
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This section is concerned with an analysis of the probability distribution 
%ax* A theoretical distribution is fitted to observed distributions 
of Rmax for various types of filtered Jimsphere data. 


Rmax '^^s calculated for each of the 150 profiles in each of the 
six data bases described in Section II. A. As illustrated in Figures 8 
through 13, the probability distribution of R„,ax accurately 

represented by a gamma distribution of the form 


Pr 



X -. 1 -,.,- 
s^r(Y) 


X 

y* x'^‘^EXP (-x/6)dx 
0 


(14) 


for 6 > 0, Y > 0 


where the parameters y and $ are functions of the mean, x, and standard 
deviation, a, of R^ax* 


Y = (x/o)^ 


(15) 


6 = ^ (16) 

The integral in Equation 14 is evaluated by a numerical approximation 
utilizing Simpson's rule. 

Table 4 contains a list of statistics pertinent to the calculation 
of the theoretical probability distributions of R^^ajc utilizing Equation 14. 
Various percentiles from the observed probability distributions of R^^ax 
during April are given in Table 5. 

As illustrated in Figures 14 and 15, y (dimensionless) and g (m/s) 
can be estimated from Xjj^(km) according to 

Y = 27.219 -4.240 In (X^ > .125 km) (17) 

g=.140X^-®^^ (18) 

where, X_ is a measure of the largest wavelength in a particular type 
of wind profile. For residual profiles (I, II and total), X^^ is 
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Figure 9. Observed and Theoretical (Gamma) Distribution of Maximum Vector 

Modulus, Rniax» Steady State Profiles During April at Cape Kennedy 
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Observed and Theoretical (Gamma) Distribution of Maximum Vector 
Modulus, Rmax> Residual I Profiles During April at Cape 
Kennedy 
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Figure 13, Observed and Theoretical (Gamma) Distribution of Maximum Vector 
Modulus, Rmax» ''n Total Residual Profiles During April at 
Cape Kennedy 



TABLE 4, mean, I, AND STANDARD DEVIATION, 0 , OF Rmax AND THE 
PARAMETERS g and y OF THE GAMMA DISTRIBUTION OF Rmax AND THE 
PARAMETER USED IN EQUATIONS 11 and 12, APRIL, CAPE KENNEDY 


Data Base 

X 

(m/s) 

(m/s) 

6 

(m/s) 

Y 

(dimensionless) 

^m 

(Km) 

Jimsphere 

42.36 

12.41 

3.64 

11.65 

48U) 

Steady 

State 

41.51 

12.67 

3.87 

10.73 

4S(1) 

Wind Bias 

38.78 

12.78 

4.21 

9.20 

43 ( 1 ) 

Total 

Residual 

7.11 

1.35 

.26 

27.91 

2.47 

Residual II 

5.77 

1.24 

.27 

21.59 

2.47 

Residual I 

2,30 

.426 

.0786 

29.30 

.42 


^^^Estimatsd; refer to succeeding paragraphs of text 


TABLE 5. PERCENTILES FROM THE OBSERVED PROBABILITY 
DISTRIBUTIONS OF Rmax DURING APRIL AT KSC 




Data Base 

50 

80 

90 

95 

99 Percentile 

Jimsphere 

41 

52 

59 

65 

76 

Steady State 

40 

52 

58 

64 

76 

Wind Bias 

37 

49 

56 

62 

74 

Total Residual 

6.9 

8.3 

8.9 

9.5 

10.3 

Residual II 

5.7 

6.8 

7.4 

7.9 

8.9 

Residual I 

2.3 

2.7 

2.9 

3.0 

3.4 
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Figure 15. Parameter B as a Function of (Cape Kennedy, April, 
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the wavelength at which the amplitude response of the high-pass filter 
is equal to .50. For wind profiles dominated by a large wavelength 
component (Jimsphere, steady state and wind bias profiles), Am cannot 
be specified from a known filter function. However, an estimate of Am 
can be obtained by assuming that the average altitude (~12 km) of Rmax 
for these profiles is Am/4; thus Am = 48 km. It must be pointed out 
that the validity of equations 16 and 17 is dependent on the validity of 
the assumption described above. Therefore, the estimation of the para- 
meters Y and 6 utilizing Equations 17 and 18 could become inaccurate for 
Aj„ greater than 2.47 km. 

A comparison of Rmax probability distributions for April and 
February for two types (TRP and II) of filtered profiles is illustrated 
in Figure 16. It is indicated that the differences between the distribu- 
tion attributable to month are small compared to differences attributable 
to wavelength range (filter type). The TRP profiles contain wavelengths 
between 90 and 420 m which are not contained in Residual II profiles. 

At the 99 percentile level, these wavelengths contribute 1.6 to 1.8 m/sec 
to the distribution of Rmax calculated from TRP data; at smaller 
percentiles the contribution is slightly smaller. 

Probability distributions of the altitude, Zmax of the maximum 
modulus calculated from April data are illustrated in Figures 17 thru 20. 
Percentiles of Zmax sre listed in Table 6. It is concluded that 

• Z is largest for residual profiles 

max 

• The distributions of Zmax for wind profiles dominated by the 
same large wavelength component (Jimsphere, steady state and 
wind bias profiles) are nearly equivalent. 

TABLE 6. PERCENTILES OF Z^ax ASSOCIATED WITH R^ax 
DURING APRIL AT .CAPE KENNEDY 




Data Base 

50 

80 

90 

95 

99 Percentile 

Jimsphere 

13.1 

13.9 

14.4 

15.2 

15.9 

Steady State 

13.1 

13.9 

14.4 

15.2 

15.8 

Wind Bias 

12.9 

13.8 

14.4 

15.1 

15.9 

Residual II 

14.4 

16.1 

16.9 

17.2 

17.7 

Residual I 

16.0 

17.0 

17.3 

17.5 

17.9 



0.06 0.1 0.2 0.S I 2 6 10 20 » 40 U 60 70 M 60 66 66 66 

6r|r,^<x} 


66,6 69.6 66.66 

tA|.M71 


Figure 16. Theoretical (Gamma) and Observed Distribution of Maximum Gust 
Modulus, R^iax* During February and April at Cape Kennedy 
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Figure 17. Observed Distribution of for Residual I and Oimsphere 
Profiles During April at Cape Kennedy 
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Figure 18. Observed Distribution of for Residual II and Jimsphere 
Profiles During April at Cape Kennedy 
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Figure 19. Observed Distribution of Steady State and Jimsphere 

Profiles During April at Cape Kennedy 
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Figure 20. Observed Distribution of Z,„ax Wind Bias and Jimsphere 
Profiles During April at Cape Kennedy 



B. Gust Vector Modulus 

A large body of statistics indicating the variation of gust 
vector modulus, R, with altitude (6 to 14 km), month (February and 
April), and filtered profile type have been generated in this study. 

An analysis of the statistics has led to the following conclusions: 

• The probability distribution of R at a particular reference 
altitude is strongly related to the strength of the large 
wavelength component in the wind profiles. A comparison 

of probability distributions of R, at KSC during April 
at 12 km, for two wavelength ranges is illustrated in 
Figure 21. A reasonable upper limit for R at 12 km in the 
wavelength range from 90 to 420 m is 3 m/sec; in comparison, 
a value of 8.5 m/sec is estimated for the wavelength range 
from 420 to 2,470m. 

t Gust modulus at 12 km in the 420 to 2,470m wavelength 

range is somewhat larger during February in comparison to 
April at KSC. As illustrated in Figure 22, the observed 
distributions diverge at the large percentiles. 

• Gust modulus calculated from total residual profiles (125 
(90 < X < 2470m) increases with altitude between 6 and 

14 km. As illustrated in Figure 23, the rate of increase 
is largest between 10 and 12 km. 

Percentiles calculated from the observed probability distribution 
of gust vector modulus from total residual profiles at various reference 
altitudes are listed in Table 7. 

C. Gust Vector Modulus Associated with Maximum Wind Speed 
and Maximum Vector Shear(ll 

Distributions of gust modulus (90 < \ < 2470m) associated with 
maximum wind modulus and maximum vector shear in unfiltered Jimsphere 
profiles are illustrated in Figures 24 and 25, respectively. Gust 
data for the distributions were obtained at the altitudes of 
maximum wind speed and at the altitudes of the top of the layer of 
maximum shear for shear layer thicknesses of 100, 500 and 1,000m. 


( 1 ) 


The equation for maximum vector shear, is 


/ r 


“max = 




max 


where Au and Av are calculated over a layer thickness, Az. 
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Table 7 


Observed Percentiles of Gust Vector Modulus for Total Residual 
Profiles (90<\%2,470m) at Cape Kennedy 


X (m/Ms) 


Pr{B<x) 


• KM 

• KM 

10 KM 

12 KM 

14 KM 

tn 

Af« 

pn 


•IV 

AIK 

•U 

A^R 

•U 

ARR 

IKBCINTI 

1 

0^ 

OJU 

aic 

137 

031 

022 

030 

047 

081 

027 


0.43 

as7 

a27 

X43 

012 

021 

047 

D.88 

OV4 

074 

$ 

au 

a7i 

aB3 

14* 

077 

0.47 

OiS 

OM 

1.02 

1.2V 

10 

0.7V 

aw 

ai2 

XTO 

1.00 

070 

1.27 

1.10 

1.47 

1.H 

20 

1tUW 

1.3V 

1.1? 

1.10 

1.40 

OM 

1.W 

1.M 

•V wn 
•mvvv 

14« 

40 

1,78 

1.M 

1,W 

1.SS 

1.V4 

1.82 

187 

140 

xw 

X71 

to 

1,17 

X04 

X10 

1.D 

117 

1.70 

X08 

IVO 

X4V 

XV8 

•0 

•V 22 

X30 

143 

tio 

lU 

1.B 

X80 

X1« 

X7V 

430 

•0 

2.IB 

2.13 

1»7 

IK 

X31 

181 

400 

4,07 

487 

448 

•0 

X12 

X40 

XM 

X13 

4.20 

xoe 

4K 

4.M 

4M 

438 

K 

1»4 

XI2 

XI2 

X42 

4.W 

XI2 

7,18 

432 

470 

7.82 

I7.E 

4.03 

4.0V 

4.2V 

423 

XV6 

4.03 

7,88 

4« 

7.13 

433 

W 

4J0 

4.38 

4.37 

4JI 

418 

4.tS 

480 

488 

470 

417 











1 

t 


Figure 22. Observed Distribution of Gust Modulus, R, Calculated from 
Residual II Profile Data (420<X<2,470m) at 12 km during 
February and April at Cape Kennedy 
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Figure 24. Observed and Theoretical (Gamma) Distribution of Gust (90<x<2,470m) 
Modulus, R, Associated with Altitude of Maximum Wind in Jimsphere 
Profiles During April at Cape Kennedy 
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Figure 25. Observed Distribution of Gust (90<X<2,470m) Modulus, R, 

Associated with Maximum Wind Shear in Jimsphere Profiles for 
Layer Thicknesses of 100, 500, and 1,000m During April at 
Cape Kennedy 
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It is indicated that the distribution of gust speed associated with 
maximum wind speed can be accurately fitted with a Gamma distribution. 
The distributions associated with maximum wind shear are approximately 
normal; the slight systemmatic upward displacement of the distributions 
illustrated in Figure 25 supports the conclusion that gusts associated 
with maximum shear tend to be slightly larger when the layer thickness 
is small. 

A comparison of percentiles of various distributions of gust 
modulus calculated in this study is presented in Table 8. 


TABLE 8, COMPARISON OF PERCENTILES FROM VARIOUS DISTRIBUTIONS 
OF GUST MODULUS CALCULATED FROM TOTAL RESIDUAL 
PROFILES (90 < X < 2,470m) DURING APRIL AT 
CAPE KENNEDY 


(m/s) 


Distribution 

Reference 

50 

80 

90 

95 

99 Percentile 

Maximum 

Fig. 13 

6.9 

8.3 

8.9 

9.5 

10.3 

Hq = 14 km 

Table 7 

4.0 

5.5 

6.4 

7.6 

9.2 

Associated with 
in Jimsphere pronie 

Fig. 24 

4.3 

5.7 

6.6 

7.5 

9.9 

Associated with 
maximum vector shear 
in Jimsphere profile 

Fig. 25 

AZ(m) 

100 

5.0 

6.6 

7.3 

8.2 

8.9 


500 

4.9 

6.3 

6.8 

8.0 

9.0 


1000 

4.6 

6.0 

6.7 

7.8 

9.0 


D. Component Gusts 

Component gusts, as defined in Section IIF, are analyzed 

here as a function of altitude (6, 8, 10, 12, and 14 km) month 
(February and April) and filter type. 

Percentiles calculated for the observed probability distri- 
butions of zonal and meridional gusts are listed in Tables 9 and 10. 
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Table 9, Observed Percentiles of Zonal Component Gust in High- 
Pass Filtered(l) Jimsphere Profiles at Cape Kennedy 


FIIRUANV (m/il 

ALTITUDS 


IKM) 

DATA 

1 

S 

10 

20 SO 

80 

90 

98 M PERCENTILE 

N 

S 

RPtI 

-iSO* 

-X19 

- 1.97* 

- 1,54 XOO 

1.49* 

1.B8 

Z27* 2.70* 

11 


TRP 

-4,10* 

-X74* 

-X31* 

- 1.90* -a4i 

1.84* 

Z32 

Z72* 130 


S 

RPII 

-X25 

- 2.03* 

-XOf 

- 1.49 X53 

1.48* 

7.90* 

ZSB* XSB* 

10 


TRP 

-X75 

-X97* 

-X40 

- 1.8B X40 

1.91* 

Z47* 

X02* AOS* 


10 

RPII 

-135* 

-XSB* 

-X17* 

- 1.82* -0.44 

1.88* 

Z18* 

3.10* A18* 

18 


TRP 

-AOS* 

-3.2B* 

-Z70* 

-Z12*-X40 

Z07* 

Z70* 

3.30* ATS* 


12 

RPII 

-5.30* 

-3 JO 

-xss 

- 1,80 aso 

ZS8* 

4,00* 

4,9S* X37* 

10 


TRP 

-6.10* 

-4.10 

-X1B 

- Z28 aoo 

X2S* 

4.80* 

XSO* 7.3S* 


14 

RPII 

-S.2S 

-X77 

-X20 

- Z24 1.2B 

Z93 

X88 

4.30 S.10 

A 


TRP 

-5.70 

-A70 

-4,07 

- 2.87 aso 

XBO 

A80 

5.30 X30 

VI 






APRIL (m/f) 





8 

RPII 

-X10 

-X25* 

-1.S3 

- 1.SB* -0.07 

1.38 

1.93* 

Z22 ZSO 

•I 


TRP 

-X75 

-X73 

-X27 

- 1.83 ao7 

1.83 

Z32 

Z70 X30 


a 

RPi: 

-X70* 

-X36 

-X13* 

- 1.89* 0.08 

1.29 

1.72 

Z14 ZS7 

n 


TRP 

-4,10* 

-X7B 

-XS2* 

- ZOB* XSO 

1.82 

1.97 

Z4S XOS 

0 

10 

RPII 

-X10 

- X18 

-1.8S 

-1.40 -X33 

1.18 

1.86 

1.9B ZSO 

A 


TRP 

-X9S 

-X78 

-Z2B 

-1.80 -X31 

1.49 

1.94 

Z30 X10 

u 

12 

RPII 

-A1B 

-xso* 

-Z92* 

-Z04* -X20 

Z37 

Z8S 

X4S ASO 

e 


TRP 

-130 

-A10 

-X40* 

-2.80* X30 

Z87 

X33 

4.03 5.30 

9 

14 

RPII 

-X50* 

-A20* 

-X37* 

-2.86* 1.S3 

XB3* 

A27* 

5.3S* ITS* 

18 


TRP 

- 7.70* 

-ABO* 

-4.12 

-X47* 1.80 

4.10* 

5.00* 

8.08* X10* 



'larger absolute value based on comparison bv month 

''UWAVELENGTH RANGE OF FILTERED PROFILES: SAi-1707 


i 

RPII 

420 - Z470 m 


TRP 

90-Z470m 
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Table 10. Observed. Percentiles of Meridional Component Gust in High^Pass 
Filtered^) Jimsphere Profiles at Cape Kennedy 


FtmUAHY InVll 

ALTITUOe 


IKMI 

DATA 

1 

5 

10 

20 SO 

80 90 

98 

M PERCENTILE 

N 

8 

RFII 

-its 

-158 

-130* 

-1.88* -130 

1.70 112 

188 

170 



TRF 

-X77 

-117* 

-170* 

-118* -140 

1.94 148 

108 

4.30 

9 

1 

RFII 

-3JS 

-183* 

-130* 

- 1.73* 0.83 

1.84* 180* 

105* 

185* 

IS 


TRF 

-4.10* 

-130* 

-173* 

-103* -115 

1.90* 180* 

130* 

4.80* 


10 

RPII 

-4.58* 

-aso* 

-181* 

-120* -140 

108* 184* 

130* 

4.90* 

18 


TRF 

-&8S* 

-4.10* 

-aso* 

-155* 100 

137* a 10* 

ITS* 

5.70* 


12 

RFII 

-8.08 

-4.78* 

-170 

- 140 180 

188* aS8* 

4.48* 

5.38* 

12 


TRF 

- 7.3S* 

-190* 

-4,70* 

- 3.00 173 

170* 4.50* 

5.30* 

8.80* 


14 

RFII 

-8.10* 

-4.70 

-190 

-3.17 -1.00 

180 180 

4,90 

6.10* 



TRF 

-iai* 


-S.2S 

-4.05 -118 

153 5.00* 

5,90* 

7.-T5 

s 





APRIL Itn/ll 





6 

RFII 

-3.50* 

-183* 

-120 

- 1.80 130 

1.83* lie* 

130* 

•4.45* 

11 


TRF 

-4.08* 

-3.10 

-157 

-1.88 -ISO 

120* 166* 

135* 

4.50* 


8 

RFII 

-ass* 

-163 

-115 

-1.58 -120 

1,55 112 

187 

117 

1 


TRF 

-aso 

-ao3 

-14S 

- 1.85 117 

1,83 155 

198 

195 

10 

RPII 

-3.30 

-185 

-125 

- 1.73 125 

1.87 120 

152 

170 



TRP 

-4.10 

-ai8 

- 173 

-100 -120 

1.M 183 

110 

4,70 

0 

12 

RFII 

-aio 

-4.S5 

-180* 

-187* -150 

1.93 3.13 

154 

178 



TRF 

-7.30 

-5.70 

-4.40 

-120* 120 

170 4.15 

4.38 

4.72 

4 

14 

RFII 

-7.70 

-as6* 

-4.91* 

-3.95* -1.00 

100* 3,80* 

4.90 

185 

10 


TRP 

-ISO 

-177* 

-6.00* 

-4.90* -120 

187* 4.70 

5.75 

13C* 



'larger absolute value based on COMf ARISON BY MONTH 
*%AVELENGTH RANGE OF FILTERED PROFILES 


RPII 420 - 2,470 m 

TRF 90 - Z470 m tAl«» 
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The distributions were calculated from RP II profiles (420 < X < 2,470m) 
and TRP profiles (90 < \ < 2,470m). The difference between the percentil 
for the two types of filtered profiles is attributable to gust in the 
wavelength range from 90 to 420 m. These differences rarely exceed 
1.5 m/s; the largest difference noted in the tabular data is 2.0 m/s 
at the 1,0 percentile of the meridional component during February, 

These differences can incre<s.,e somewhat at extreme percentiles (<1, >99) 
not given in the tables. As indicated in Table 11, the largest observed 
extremal in the go to 420m wavelength range at the same reference 
altitudes is -2.61 m/s (April meridional component at 12 km). In 
comparison, the largest extremal in RP II and TRP profiles is 9.96 
and 11.70 m/s, respectively. 


The asterisks in Tables 9 and 10 indicate the larger absolute^ 
value based on a February to April comparison for equivalent percentile, 
altitude and filter type; the 50 percentile was not included in the 
comparison. The number of percentiles, N, for a particular altitude 
and month that exceed the corresponding percentiles for the other 
month is a measure of the relative strength of gusts for the two months. 
The following criteria were used to evaluate the relative strength of 
gusts: 


Relative 


N 

Strength 

Tendency 

0, 1, 2, 3 

Smaller 

Strong 

4. 5, 6 

Smaller 

Moderate 

7, 8, 9 

Equal 

None 

10, 11, 12 

Larger 

Moderate 

13, 14, 15, 16 

Larger 

Strong 

February gusts have a moderate to strong tendency to be larger at 6, 8, 
10 and 12 km for the zonal component and at 8, 10 and 12 km for the 
meridional component. April zonal gusts have a '.trong tendency to be 
larger at 14 km; April meridional gusts have a moderate tendency to be 
larger at 6 and 14 km. 


The probability distributions of the gust components are "S" 
shaped and cannot be fitted to a high degree of iiccuracy with a normal 
distribution. The probability distributions of the gust components 
during April at 12 km for RPI (90 < X < 420m) profiles (Figure 26) 
and F(P II (420 < X < 2470m) profiles (Figure 27) illustrate the "S" 
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shape. Distributions of this type have a probability density function 
which has two maxima. If the maxima are nearly symmetrical about zero, 
the density function of the absolute values has a single maximum and 
the shape of the probability distribution can be accurately approximated 
by a ganma distribution. The form of the gamma distribution for gust 
components is the same as Equation 14 (Section 1 1 1. A) with either |ul 
or lv| substituted for 

An example of the observed and theoretical (gamma) distribution 
of lu] during February and April at 12 km is illustrated in Figure 28. 

A list of statistics pertinent to the calculation of theoretical 
distributions for juj and [vj for TRP profiles is given in Table 12. 

2. Component Gusts Associated with Maximum Wind Speed and 
Maximum Wind Shear 

Component gusts associated with maximum wind speed and 
maximum vector wind shear in unfiltered Jimsphere profiles were derived 
from April data according to the procedure outlined below: 

e The altitude of maximum wind speed was identified. 

• The altitude of the top of the layer containing the maximum 

vector shear was identified for shear layer thicknesses 
of 100, 500, and 1,000m. 

§ The altitudes described above were used for determination 
of gust component amplitudes in high-pass filtered (TRP) 
profiles (90 <X< 2,470) according to the procedure described 
in Section II, A. 5. This yielded four sets of component 
gust data (three associated with maximum shear and one 
associated with maximum wind speed). 

The distributions of gust components associated with maximum shear are 
illustrated in Figures 29 (zonal) and 30 (meridional). All the distri- 
bution have an "S" shape which is characteristic of a probability density 
function with two maxima. It is indicated that 70 percent of the zonal 
component gusts associated with the larger vector wind shears (aZ= 500 
and 1,000m) are negative. Meridional component gusts do not exhibit 
such a strong bias. 

Distributions of gust components associated with maximum wind 
speed are illustrated in Figure 31. These distributions also exhibit 
and "S" shape. A large percentage (>94^) of zonal component gusts 
associated with maximum wind speed are positive; a smaller percentage 
(61%) of meridional gusts are also biased but in the opposite sense 
(negative). Therefore, it could be concluded that the contribution of 
gust in the wavelength range from 90 to 2,470m to the maximum wind 
speed in a Jimsphere profile is more likely to be positive than negative. 
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Table 12. Mean, Standard Deviation and Parameters y and 6 of the Gamma 
Distribution for Absolute Component Gust in Total Residual 
Profiles (90<X<2,470m) at Cape Kennedy 


FEBRUARY 



6 KM 

8 KM 

10 KM 

12 KM 


14 KM 


|U| 

M 

luj 

|v| 

iul 

|v 1 

Iul 

Iv| 

Iu| 

M 

MEAN (m/i) 

1.62 

1.86 

1.72 

1.81 

1.90 

2.25 

263 

3.07 

288 

3.57 

S. D. (m/i) 

0.90 

0.97 

0.94 

1.07 

1.08 

1.33 

1.72 

1.80 

1.58 

200 

P (m/i) 

aso 

0.51 

a512 

j 

0.63 

0.61 

0.79 

1.13 

1.06 

0.87 

1.12 

7 

3.24 

3.68 

3.35 

2.86 

i 

3.10 

2.86 

234 

2.91 

1 

3.31 

3.19 





APRIL 






MEAN (m/s) 

1.S9 

1.87 

j 

1.59 

1.72 

1.45 

1.81 

2 37 

273 

237 

3.69 

S. D. (m/s) 

0.86 

1.01 

0.84 

0.99 

0.82 

1.10 

1.32 

1.60 

1.90 

2.07 

0 (m/s) 

0.47 

0.55 

0.44 

0.57 

0.46 

0.67 

0.74 

0.94 

1.07 

1.16 

7 

3.42 
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3.58 
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222 
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3.15 
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Figure 28. Observed and Theoretical (Gamma) Distribution of Absolute Value 
of Zonal Component Gust (90<X<2,470m) at 12 km During February 
and April at Cape Kennedy 
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Figure 29. Distribution of Zonal Component Gust (90<X<2,470m) Associated 
with Maximum Vector Shear in Jimsphere Profiles for Layer 
Thicknesses, AZ, of 100, 500, and 1,000m During April at 
Cape Kennedy 
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Figure 30. Distribution of Meridional Component Gust (90<X<2,470m) Associated 
with Maximum Vector Shear in Oimsphere Profiles forLayer 
Thicknesses, AZ, of 100, 500, and 1,000m During April at Cape 
Kennedy 
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Figure 31. Distribution of Component Gust (90<X<2,470m) Associated with 
Maximum Wind Speed in Jimsphere Profiles During April at 
Cape Kennedy 
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Extreme percentiles of component gust (90 <X< 2,470m) distribution 
are compared in Table 13. A strong similarity between all but one of 
the distributions is noted; the exception is the relatively small negative 
zonal component gusts associated with maximum wind speed. Gusts associated 
with maximum wind shear or maximum wind speed are not significantly 
larger than those at a reference height of 14 km. 

F. Gust Length 

Gur>t length. L, as defined in Section II. A. 5 is analyzed 
herein as a function of altitude (6, 8, 10, 12, and 14 km) and month 
February and April) for TRP data (90 <X< 2,470m). 

Percentiles calculated from the observed probability distributions 
of L for zonal and meridional component gust are listed in Table 14. It 
is indicated that thore are no significant trends in the percentiles as 
a function of altitude, month and component, The observed and theoretical 
(gamma) distribution of L at 12 km is illustrated in Figure 32., The 
parameters for the gamma distribution illustrated in Figure 32 were 
calculated from meridional component gust length data at 12 km during 
April. The gamma distribution can provide an accurate estimate of the 
observed distribution. The observed distribution is bounded at the large 
gust length end of the distribution by the filtering process used to 
derive TRP profiles; therefore, the gamma distribution, which is un- 
bounded at large gust length, is not valid outside the gust length range 
of the observed distribution. From this it is suggested that a beta 
distribution may provide a more accurate approximation. 

The beta distribution is of the form 


I 

Pr IL < Tj = f 


(a+6+1), 
aJg; 

where t = x/x^g^ (^Itll) 


(19) 

( 20 ) 


The parameters a and 3 are estimated from the sample mean, X, 
and standard deviation, a, according to 


a = Kb - 1 

(21) 

* 

S = b - 1 

(22) 
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Table 13. Observed Extreme Percentiles of Component Gust (90<X<2, 470m) 
Distributions During April at Cape Kennedy 


u (m/i) 


DISTRIBUTION 

AZ (m) 

1 

5 

10 

90 

95 

99 

PERCENTILE 

ASSOCIATED WITH 

lOO 

-7.12 

-6.40 

-121 

6.67 

7.70 

136 


MAXIMUM VECTOR 
SHEAR IN JIMSPHERE 
PROFILE 

500 

-8.06 

-6.20 

-5.69 

126 

199 

7.86 



1000 

-7,12 

-6.11 

-141 

111 

181 

7.97 


ASSOCIATED WITH 
MAXIMUM WIND SPEED 
IN JIMSPHERE PROFILE 


-170 

-0.30 

1.60 

110 

170 

110 


REFERENCE HEIGHT 
- 14 KM 


-7.70 

-4.90 

-4.12 

100 

105 

110 





w (m/f) 





DISTRIBUTION 

AZ (m) 

1 

5 

10 

90 

95 

99 

PERCENTILE 

ASSOCIATED WITH 
MAXIMUM VECTOR 

100 

-117 

-7.64 

-185 

111 

196 

113 


SHEAR IN JIMSPHERE 

500 

-9.09 

-7.38 

-6.19 

119 

7.38 

137 


PROFILE 

1000 

-109 

- 7.22 

-171 

172 

7.71 

172 


ASSOCIATED WITH 
MAXIMUM WIND SPEED 
IN JIMSPHERE 
PROFILE 


-150 

-165 

-180 

4.30 

155 

7.10 


REFERENCE HEIGHT 


-150 

-177 

-100 

4,70 

175 

130 
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Table 14* Observed Percentiles of Gust Length, L, in TRP (90<X<2,470m) 
Profiles at Cape Kennedy 


L (m), APRIL 



6 KM 

8 KM 

10 KM 

12 KM 

14 KM 

PERCENTILE 

u 

V 

u 

V 

u 

y 

u 

V 

u 

y 

1 

38 

25 

55 

13 

2S 

19 

25 

25 

15 

9 

5 

90 


138 


107 


88 


81 




131 


125 


79 


105 

^ * 

42 

10 

150 

188 

194 

225 

141 

163 

167 

180 

150 

100 

20 

300 

313 

350 

365 

238 

309 

350 

310 

380 

270 

50 

609 

589 

639 

694 

580 

650 

625 

610 

660 

556 

30 

967 

SS6 

870 

950 

950 

900 

906 

875 

879 

825 

90 

1025 

1027 

1088 

1213 

1117 

1017 

1067 

1075 

984 

9;’.5 

95 

1142 

1144 

1235 

1363 

1275 

1175 

1159 

1213 

1131 

996 

99 

1588 

1388 

1235 

17S 

1725 

1375 

1663 

1738 

1525 

1425 





L (ml. FEBRUARY 







6 KM 

8 KM 

10 KM 

12 KM 

14 KM 


PERCENTILE 

u 

V 

u 

y 

u 

y 

u 

V 

u 

V 

1 

54 

25 

25 

10 

37 

19 

13 

9 

10 

10 

5 

97 

135 

88 

47 

96 

94 

63 

42 

47 

47 

10 

167 

250 

134 

94 

139 

170 

167 

78 

117 

100 

20 

317 

405 

300 

317 

292 

350 

250 

189 

315 

263 

50 

534 

659 

360 

618 

613 

644 

560 

593 

636 

584 

SO 

800 

915 

880 

890 

928 

960 

838 

979 

959 

916 

90 

1063 

1110 

1050 

1050 

1134 

1200 

ir>oo 

1138 

1086 

1050 

95 

1213 

1215 

1275 

1175 

1225 

1413 

1257 

1263 

1231 

1155 

99 

1563 

1525 

1575 

1588 

1563 

1583 

1525 

1613 

1725 

1375 


SAI-9U6 
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Observed and Theoretical (Gamma and Beta) Distributions of Gust 
Length Calculated from TRP Profile Data (90<X<^2,470m) at 
Cape Kennedy 
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(23) 


K = 


Vl 


max 


1-X/X, 


max 


K-(g/X (Kn)^ 


(24) 


The beta distribution illustrated in Figure 32 was also calcu- 
lated for the meridional component gust length data at 12 km during 
April. The value of X^ax (3,125m) chosen is one-half the wavelength 
at which the high-pass filter amplitude response for TRP profiles is .10. 

As illustrated in Figure 32, the beta distribution does not seem to fit 
the observed distribution with greater accuracy than the gamma distribution. 


The parameters of the gamma and beta distributions of gust length 
are listed in Table 15. 


S. Prediction of Monthly Wind Statistics 


The year=to-y8ar variability of monthly winds aloft statistics 
at Cape Kennedy has been documented in a previous study (Ref. 8). 

The bivariate normal statistics during April at Cape Kennedy during 1958 
and 1967 listed below illustrate the magnitude of the variability. 


(m/s) (m/s) (m/s) (m/s) 

1958 45.89 19.78 -6.21 16.70 .0176 

1967 14.14 12.75 -9.38 7.90 -.1610 


For the purpose of launch planning, it is desirable to optimize the 
selection of monthly wind statistics so that conditions associated with a 
particular launch are more accurately represented. The selection could be 
based on 30-day forecasts of synoptic weather patterns issued by NOAA. 

An examination of the synoptic patterns during 1958 and 1967 
indicates that there, is a correlation with certain wind statistics. 

As indicated in Table 16, the persistant high-pressure ridge over the 
southeastern U.S. during 1967 resulted in unusually warm and dry conditions', 
in contrast the synoptic pattern during 1958 favored above normal 
precipitation and near normal temperature. The rank of each variable 
during the fifteen-year period (1956-70) is included in the table; for 
example, the average precipitation of .08 inches in 1967, the smallest 
in the 15-year period (1956-70), has a rank of 15. The data in the table 
indicate that the wind statistics are positively correlated with 
temperature; the correlation with precipitation seems to be higher. 
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The foregoing examination is the basis for an attempt to predict 
total wind variability during a future 30-day period at Cape Kennedy 
given the sta^idard NOAA 30-day forecast (updated every 15 days) of 
temperature and rainfall. The method is based on the hypothesis that 
total wind variability at Cape Kennedy is correlated with mean temperature 
and rainfall over a climatic division (area) near Cape Kennedy. The 
total wind variability, Oj, is defined by 

“T 0 / + (25) 


where, o' and are the standard deviations of the zonal and 
meridional wind'^components during a monthly reference period. 


The standard NOAA 30-day forecasts of temperature and percipitation 
are of a qualitative nature and cover large areas. Precipitation predictions 
are in three broad categories (light, moderate and heavy) each representing 
one-third of the monthly values that have been observed in the past; 
similarly, there are five categories (much below, below, near normal, above 
and much above) for temperature; one-eighth of the observed temperatures 
that fall in each of the much above and much below normal categories, one- 
quarter in each of the above and below normal, and near normal categories. 

The class limits have been determined by NOAA for each month and station 
from available uniform records. The class limits for Orlando, Florida, 
during April are given below along with values that have been chosen to 
represent each category, 

_ Representative 

Average Temperature (T) Class Limits (°F) Values (°F) 


Much Above 
Above 

Near Normal 
Below 

Much Below 
Precipitation (P) 


73.4 < T 

71.9 < T < 73.4 

70.5 < T < 71.9 
69.0 < T < 70.5 

T < 69.0 


Class Limits (in) 


73.4 

72.7 
71.2 

69.7 
69.0 

Representative 
Values (in) 


Heavy 

Moderate 

Light 


4.18 P 4.18 

1.52 D < 4.18 2.85 

P < 1.52 .76 


The total wind variabi_lity, a, is estimated by substitution of the 
representative values of T and P into the relation 

a = 78.68 + 1.50P - .88T (26) 
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The coefficient of multiple correlation between a, P and T is .80. 
Equation (26) was calculated by the least squares method from data for the 
period 1956-70; Cape Kennedy wind data and mean values of P and T obtained 
for the North Central Climatic Division of Florida were used. 

The correlation between the observed and calculated a at 12 km 
during April utilizing P and T observed during the period 1956 thru 
1970 is illustrated in Figure 33. It is shown that even if P and T 
are known, or predicted exactly, significant differences exist between 
the calculated and observed wind variability. These results illustrate 
the best that can be achieved with this technique. Since the NOAA 30-day 
forecasts of T and T are not perfect, a further degradation of the 
correlation between observed and predicted a is expected when the fore- 
casted values of F and T are used in Equation (26); this is illustrated 
in Figure 34. It is evident from these results that the NOAA 30-day 
forecast is not sufficiently accurate to be utilized in Equation (26) 
for the prediction of o on a long term basis. 
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(S/Ul) o 


(+) OBSERVED AT 12 KM OVER CAPE KENNEDY DURING APRIL 

(O) CALCULATED FROM EQUATION (26) UTILIZING OBSERVED P AND T. 



YEAR SAt-5690 


Figure 33. Observed and Calculated a 12 jcm During April at Cape 
Kennedy Based on Observed P and T 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


Probability distributions of a number of characteristics of detailed 
wind profiles measured at Cape Kennedy, Florida, have been shown to be 
accurately represented by a gamma distribution. These characteristics include: 

• Maximum vector modulus of filtered and unfiltered Jimsphere 
wind profiles 

t Gust vector modulus at a reference height 

• Gust vector modulus associated with maximum wind speed in the 
unfiltered Jimsphere profile 

• Zonal and meridional component gust at a reference height 

• Gust length 

There is an apparent inconsistency in the conclusions concerning 
component gust and gust vector modulus at a reference height. If the 
component gusts are bivariate gamma, then the gust vector modulus would not 
be univariate gamma as indicated here. Lacking a complete knowledge of the 
probability distribution of the modulus of a bivariate gamma distributed 
vector, it can only be concluded that the univariate gamma is a good approximation 
of that distribution for this data set. The validity of the univariate 
approximation may be proven if a better understanding of the bivariate gamma 
distribution and the required simplying assumptions can be developed. It is 
recommended that methods for calculating the probability distribution of the 
modulus of the bivariate gamma distributuon be devised. 

The results of the study provide a basis for the further development 
of idealized models for use in flight simulations. These idealized models 
by themselves cannot represent real wind profiles which have a significant 
random component. It is recommended that a model consisting of the attributes 
of both idealized and random gust models be developed. 
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